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1. INTRODUCTION 

How the chaotic matter did get organized into compact astrophysical objects, such as stars, when the Universe was 
created? How did the Sun's rotation get slowed down (a central problem of cosmogony)? What is the physics behind 
quasars' shining? The answers of all those unresolved problems of contemporary physics go back to the problem 
of the effective viscosity of weakly magnetized plasmas in shear rotating flows. For half a century we have faced a 
kinetic problem — how to calculate an effective viscosity — a problem that is at the core of the machine for making 
stars. This longstanding problem has already been approached in so many ways — any proposal-writing astrophysicist 
has already published his/her view and the literature of analytical works and numerical simulations is overwhelming. 
Yet the problem is still unsolved and has not lost its attractiveness. Here we will give it a try, too. 

2. MODEL AND MHD EQUATIONS 

Our starting point are the conservation laws for energy and momentum for an incompressible fluid with mass density p 

d,{pVi) + d,{Uik)^0, (1) 
-{li—+p£+—]+dwq^O, (2) 



dt \ 2 2^0, 
divV = 0, p = const, (3) 

where we have for total stress tensor n and heat flux q respectively 

n,, - pViVk + P5,k -'l[^^+^)-J[-^ {B,Bk - ^5ikB^^ , (4) 

— V-V-C7'-K-Vr+— [Bx (VxB)]-^^(BxrotB), (5) 
2 / IM) Mo 
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where e is the internal energy per unit mass, w is the enthalpy per unit mass, (7^ = 77 (5,Vjt + d^Vi) the viscous part of 
the stress tensor for an incompressible fluid, 77 is the viscosity, K is the heat conductivity, T is the temperature, p is 
the Ohmic resistivity. The formulas are written in SI, for a transition to Gaussian system we substitute /io = 4-n and 
£0 = 1 /4-n, i.e., expressions are written in an invariant form. 
For the magnetic field's energy density rate of change we have 

-LIb' = 1b.^ = 1b.[Vx(VxB)-v„,Vx(VxB)] 

= J_V.(Bx (VxB)-VmBx (VxB)). (6) 
Mo 

We calculate the divergence of the total stress tensor Eq. (4), and using that 

dk— (BiB, - 1^25 ^ = 1 ("b . VB - i VB^) = — (B X rotB) = j x B, (7) 
Mo V 2 y Mo V 2/^0 

we obtain the equation of motion for an incompressible plasma, 

p<9,V=-V-VV-VP+jxB + 77V2v. (8) 

To close the set of equations we need to use Ampere's law, V x B = Moj, and generalized Ohm's law, E + V x B = Tjj, 
and supplement them with Faraday's law 

5B 

— = -VxE. (9) 
at 

For the second set of MHD equations we have 

= Vx ( VxB-— VxB). (10) 



dt \ Mo 

The MHD equations for an incompressible fluid p = const, in homogeneous magnetic field Bo, shear flow with rate A, 
and angular velocity SI =A(oe^- = Afi) are 



VxB 



pD,V = -VP+ |^j= — — j xB-2pQx V + pVkV^V, (11) 
D,B = B-VV + VmV-B, divV = 0, divB = 0, (12) 

where D, = (9, + V • V is the substantial (convective) derivative, P is the pressure, j is the current density and Vk is the 
kinematic viscosity. The magnetic diffusivity Vm = eoc^P is expressed by the Ohmic resistance p and £0 = 1 /Moc^- In 
order to obtain a linear system of dimensionless MHD equations we use the following Ansatz for the velocity V, the 
magnetic field B, the wave vector Q, and the pressure P 

V(r,r) = V.,hear(r)+Vwave(f,r), Vwave = iVA^VQ(T)e''3-^, Vshear=Axev, (13) 

Q 

B(f,r) = Bo+Bwave(f,r), Bwave(r,r) = BoEbQlTje'^^-^, (14) 

Q 

P(f,r) = Po + ^'wavel^r), Pwave(f,r)=pyi£^Q(T)e'Q-^. (15) 

Q 

where the sums are actually integrals with respect to 3-dimensional Eulerian wave-vector space with independent 
coordinates g., , Qy, Qz 

AQAQyAQ, f,ifQ\ f d^e 



Q 



i.e., the sum is a short notation for Fourier integration with omitted differentials, integral limits and In multipliers. 



For the static magnetic Bq field with magnitude Bq = ^IBq^,+Bq_ we suppose a vertical Bq- and an azimuthal B^y 



components parameterized by an angle 9 and an unit vector a. We also assume that the Alfven velocity Va is much 
smaller than the sound speed Cg 



Bo = BqU, a = (0, a,, = sin 0, = cos 0), V/ 



Bo 



Bp 

/Mop' 



(17) 



In the equations above we used the dimensionless space-vector 

X 

I r Mr 

A = 




Va 




(18) 



and its dimensionless wave-vector counterpart Q. Here, A is the characteristic length of the system which we suppose 
to be much smaller than space inhomogeneities, e.g., the accretion disk thickness. 



3. WAVE- VECTOR REPRESENTATION 
3.1. Linear terms 

We have a space homogeneous physical system and indispensably its modes bear the character of plane waves. The 
purpose of the present section is to find the Fourier transformation of the all the terms in the MHD equations (13), (14) 
and (15). 

Let us start, for example, with the pressure. According to Eq. (15) we have 

y = ? / ^Q^")*^^'"'' (0 = iAyAi:QPQ(T)e'Q^ (19) 
VP t' VP 

— ) = J c-'^-^—d^X =iAVAQPQiT). (20) 

Analogously, according to Eq. (13) and Eq. (14), for the partial time derivatives we obtain 

o',V = iAyA5^e'Q-^5,VQ(T), .#(5,V) = ;-AyA^,VQ(T), (21) 
Q 

a,B=ABo£e'<3-^a,bQ(T), #(5,B)=ABo5,bQ(T). (22) 
Q 

More complicated are the Fourier transformations of the expressions, containing the shear flow Vshear = Va^Cv and 
the wave variables 



Vshear -VVwave = -AVa J ^^-^XQyyiT) 
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(23) 



Vshear -VBwave = iABo / e'^^'^XOM^)^^ ■ (24) 

Let variables Vwave or Bwave be presented by their Fourier components \ff{X) = T.q^^^^Wq ^^'^ Vq = ■^iwOQ)- Our 
task is to derive the Fourier transformation ^(X\//(X)). Using that Xe'*-*'^ = — i^ge''^ '^ and the Gaussian theorem, 
Jyd^QdQ ~ fgf^dS applied for the whole volume in wave-vector space and its boundary dY we can make the 
partial integration 

Xv/(X) = X^e'Q-^rQ = -i^ {(5Qe'Q-^)v/Q = -e'^'^^QV/Q + [cQ^^^qV/qI | = ^e'^'^i^QV/Q, (25) 

Q Q '- ^ Q 

in the Umit 

lim (eVQ)=0. 



In such a way we derived the well-known in quantum mechanics operator representation X = idQ and derived the 
Fourier transformation 

#(XvA(X))=iaQV/Q. (26) 
This expression is analogous to the Fourier transformation of the V-operator 

#(Vx) = iQ, (27) 

and yields 

^{XtyVx)^-QydQ^. #(Vshe»-V) = -Aej%,., Vshear - VaXCj. (28) 

Those relations give that 



D;'^'=^^ = d,+ V,hear ' V = ^, +AXd,^ = A {of = d, - Q,dQ^ ^ d, + Ushear(Q) ' ^q} • (29) 



In other words Fourier transformation of a linearized substantial derivative is again a linearized substantial derivative, 
but only in the wave-vector space. For this purpose we introduced the field of shear flow in the wave-vector space 
Usheai (Q) = —Qy^x', Confer this result with Vshear/^A — Xty Returning back to the velocity and magnetic field we 
arrive at 

^[(5r+V,hear-V)Vwave] =iAyA[5.+U,hear(Q)-5Q]VQ, (30) 
Vshear-V)Bwave] = ABo[<5t + Ushear(Q) ' ^qlbq . (31) 

For the derivation of these equations we used Eq. (28) and according to Eqs. (21) and (22) J^{di) = Ad^. 

Very simple is the Fourier transformation of the dissipative terms which is reduced to the properties of the Laplacian 

v,V^V = -v,.y. / VQ(r)G^e'Q-|| - v,/ VQ(r)G^e'<^-|| = -iAy.v( / G^VQ(r)e'Q-|| , 



#(VV2V) - -vLe'v(T), — #(VmV2B) = -V;e2b(^), V^^-^Vk, vL^7j2^m. (32) 

iAVa ABq Vt Vi 



Hereafter for all terms coming from the velocity equation (11) we will separate a factor IAVa and for all terms from 
Eq. (12) we will separate a factor ABq - Those factors will be common for the final equations in the wave-vector space. 
For Coriolis force density per unit mass — 2£2 x V, we derive 



Ax + iyA/v,Q(T)eiQ-x||, 



-2ax V= -2A»(-y,e, + y,ey) =2AC(J | -V, j = 2A» |^ ^iVf,,\ v,^q{x)^^-^^ 

#(-2£l X V) = iAyA2a)(v,,Q(T)e, - v,,Q(T)e,). (33) 

The centrifugal term 2(oA^x is irrelevant for the wave amplitude equations. 

Furthermore, we calculate Lorentz force per unit mass j x B, and using Bq/Mo = P^a have 

'^^)"t) - ^/(Qx^,x(.„,,e,+«.,e,e.»''||.MV./[(Q.b,W)xa,0, 

# ( ( y^^\ = IAVa [(Q x bQ(T)) x a] . (34) 



Mo y P 

We have also other two zero terms having no influence on the wave dynamics. From momentum equation (11) and 
the magnetic field equation (12) we have 

Vshear ' VVshear = Axty ■ VAxCv = A^Xty ■ t^ty = , (35) 

Bo • VVshear = ABo(«,.ej, + a,e,) • e.,ev = . (36) 



For the last linear terms we have 



Vwave • VVsheai- = l4Va j Vq(t) • ^^^y^^'^ -^§^ ^ ^(Vwave ' VVshear) = iAyAV;,,Q(T)ev , 



(37) 



Bwave-VVshe^ ^ ABq J b,,Q{x)eye'^-^ , #(Bwave- VVshear) =ABo^,-,Q(T)ev, (38) 



Bo-VVwave = /(a-iQ)VQ(T)e'Q-^y:^, .#(Bo-VVwave) = -AZ?o(a-Q)VQ(T). (39) 



^VaBo , ^ lOX 

(2;r)3 

All linear terms are well-known from previous investigations of MHD waves in magnetized shear flows. In the next 
subsection we will derive the nonlinear terms describing the wave-wave interaction coming from the convective time 
derivative. 



3.2. Nonlinear wave-wave interaction 

In order to derive the nonlinear term in the momentum equation (11) we calculate Vwave ■ V Vwave using 
l/A 

Vwave • VVwave = ^ iVAv(T, QOe'^' ^^ • V^iVAvCr.Q'Oe''^"-^ = -iAVAli: VQ, • Q'Ve'(Q'+Q") '' . (40) 

Q' Q" Q' Q" 

For the sake of brevity in the last terms we will omit the time argument T and write the wave-vector argument Q as 
index. We make a Fourier transformation and obtain 

# (Vwave • VVwave) = J d^X^-'^-^ (Vwave ' VVwave) (41) 

= -iAVAEE^Q' -Q'V^ / Q +Q -Q \ ^ _L4yA^VQ_Q,VQ, Q. 

Q' Q" \ / Q' 

The velocity and magnetic fields 

Vwave(T,X)=iyA/eQ'^v(T,Q)^, Bwave (t, X) = Bo | cQ'^Mt, Q) ^ , (42) 

have to be real, hence the Fourier components should be odd for the velocity and even for the magnetic field 

V Q= vq, b Q=bQ. (43) 
Analogously, for the Fourier component of the wave-wave interaction part of the Lorentz force jwave x Bwave we obtain 

# ( — (V X Bwave) X Bwave ) = / d^Xe''^^ f — (V X Bwave) X Bwave ) = iAVA^ (Q' X bg,) X bg^Q,. (44) 

VMoP / J VMoP / Qi 

In such a way we derive the Fourier component of the nonlinear term of the momentum equation 

NvQ ~ ■ .1, ^ ^ Vwave ' VVwave H (V X Bwave) X Bwave ) (45) 

iAVa V MoP / 

= L [>'Q-Q''*'Q' • Q + (Q' X bgO ^ Q'] ■ (46) 
Q' 



Similarly, for the other nonhnear terms Bwave ■ VVwave and Vwave ■ VBwave we have 

#(Bwave-VVwave) = ^ d^Xe-'^'^ (Bwave ' VVwave) (47) 

Q' Q 

# (Vwave -VBwave) = / d^Xe-'^'^ (Vwave • VBwave) (48) 



Q' Q" \ ^■"^ / Q> 



Those terms participate in the equation for the magnetic field. For their difference we have 

N/7 Q — , ^ (Bwave ' V Vwave Vwave ' VBwave ) (49) 

ABq 

= E [^Q' (Q Q') bg Q' bQ, (Q Q') vq q ] 

Q' 

= -QX^(Vq/ XbQ_Q/) 

Q' 

As the function in r-space 

rot (Vwave X Bwave) ~ I^wave ' VVwave Vwave ' VBwave (50) 

has zero divergence 

div [rot (Vwave X Bwave ) ] = (51) 

its Fourier transform is transversal Q • N^.q = and automatically N^q = N^ q. 

In order to merge the so derived nonlinear terms in a next subsection we will rederive the linear terms in the 
Lagrangian wave-vector space. 

4. ELIMINATION OF PRESSURE IN THE FINAL MHD EQUATIONS 

It is common in MHD to formally seek the limit of a particular expression for infinite sound speed Cg — > °°. Due to the 
complexity of the problem this standard approach for consideration weak magnetic fields when Va <C Cj is inapplicable 
to our problem and we have to look for direct elimination of the pressure. After substituting Fourier transformations 
in Eqs. (11) and (12) we obtain 

((9T + Ushear-<5Q)VQ -V^,Qev + QPq + [(Q X bg) X tt] + 2a)(v,,Qe.v - V^^qC,.) - v^g^ (52) 

+L [vq-q'Vq' ■ Q + (Q' X bQ/) X bQ_Q/] , 
Q' 

(5T + Ushear-5Q)bQ b^.Q^,. - {Q- a)\Q - V^Q^hQ - Q X Y^{\q, X hQ_Q,) . (53) 

Q' 

For the sake of brevity we introduce 

= Qy ^ - e,e, • vq + [(Q x b) x a] + 2(0 x vq + vIQ^\q + ^ [vq_q,vq/ • Q + (Q' x bg/) x hQ q,] . (54) 

Ofjx Q/ 

Then the equation for the velocity can be rewritten as 

^tVq=PqQ + ^q. (55) 
In order to express the pressure, we multiply both sides of this equation by Q 

^T(QvQ) = e2pQ + Q.^Q. (56) 



The incompressibility condition Q • vq = gives for the pressure the solution of the Poisson equation 
^__Q^ = --i2{2Q-e,e,-VQ + 2a)xvQ + Q-[(QxbQ)xa]} 



1 



2 L {Q ■ '*'Q-Q'Vq' ■ Q + [(Q' X bgO X bQ_Q,] • Q} , (57) 
Q' 

where we used the obvious vector relations 
Q-eve;,-VQ = 2y,-ej., Q • (O) x vq) = ©(GvV.v - G.-Vv), (Q x bg) x a = (Q • a)(Q -b) - 0^(5 • a) . (58) 
This formula for the pressure we substitute in Eq. (55) which takes the form 

^TVQ-^^=^Q-^^^Q=n^'2^Q, (59) 

where 

n^Q^ll_Q^ = ll_„«„, (60) 

G^ G 

is the projection operator which applies to the part of a vector, perpendicular to the wave-vector In other words, the 
elimination of the pressure conserves the perpendicular part of the Fourier component of the force in the used 
dimensionless variables. Equation (59) means that the velocity field remains orthogonal to the wave vector If in the 
beginning Q • vq(to) = 0, the evolution gives that Q • vq(t) = for every T > Tq. 
Using that for the velocity as apphcable for every orthogonal vector 

n^a,VQ = a,VQ, n% = VQ (61) 

we can rewrite Eq. (59) as 

n^(5TVQ-^Q)=0. (62) 
In order to take into account the Qydy^^ /dQx term in Eq. (54) we use the obvious relations 

Now we represent the projection of the advective term Ushear • <5qVq as 

n^Q(Ushear ' <3qVq) = ~QyJ^ " «vnv,,,Q (64) 

and arrive at the momentum equation in the form where the projection operator exists explicitly only in the nonlinear 
term 

{d-c +Ushear ' <5q) vq(t) = -Vr.QCy + InyHv^^Q + 2c(jn(nja'.,. q - «xVv,q) + 2(0 X VQ + (a • Q)bQ 

-V^G'vq +n^Q^ [vq_Q'Vq, • Q + (Q' X bgO X bQ_Q,] , (65) 

Q' 

(<3t + Ushear ' ^q) bQ(T) = b^^qty - (Q • a)VQ - Vj^^hf^ - Q X ^(vq/ X bg^Q/), (66) 

Q' 

Vq(To) = n^VQ(To), bQ(To) = n^bQ(To). (67) 

For numerical calculations the incompressibility conditions n • bg = and n • vq = can be used as a criterion for the 
error. 

Using the relation 

[Ushear ' ^qbg (t)] • Q = Gv^-.Q, (68) 



one can easily check that the equation for the evolution of the magnetic field, Eq. (66), can also be presented as the 
evolution of its part, perpendicular to the wave-vector 



<9i:bQ + n^Ushear ' <5QbQ = n^CyC, • bg - (Q • a)\Q - v'^Q^bq - Q X ^(vq/ X bQ_Q/). 

Q' 

Together with vq = H^vq this equation automatically gives bg ~ n^bg and divB = 0. 
In the matrix form the set of MHD equations reads as 



(69) 



where 



and 



M : 











QxlQ'(VQXbQ Q/) 

[vQ_Q'VQ/-Q + (Q'xbQ/) 


X bg Q 


1 






-vLQ^ 







-<Q^ 


~Qa 








-Qa 






-Qa 


Qa 






Qa 






Qa 


lnynx[(£) + \)- 
2nyny{a)+ 1) — 
2;iv;i;(£(J + 


-KQ' 

(2(0 + 1) 
1) 


-2nxnx(0 + 20} 
-2nxny(0 - v[Q^ 
—2nxn,(0 







with Qa=Q a. 
The matrix can also be represented as 



M 



n.v 








riy 








«, 









Mhb I Mfc, 
Kb I M„v 



1 




-Qat M 



bb 







{nl + nl) 


-rixriy 




-rixn^ 


/ 




1 




V 





(70) 
(71) 





( bx 


\ 




by 






bz 






Vx 
V,, 








/ 



(72) 



(73) 



(74) 



(75) 



With the help of the matrices in this representation, in the next section we will make Lyapunov analysis of the linearized 
MHD equations. 



5. LYAPUNOV ANALYSIS OF THE LINEARIZED SET OF EQUATIONS IN 

LAGRANGIAN VARIABLES 

For small Qy we may neglect the advective term Ushear • <5q = —Qydg^ - Then the linearized MHD equations take the 
form 

D^^=M'P, (Q,Q)-^ = 0. (76) 

To perform an instability analysis we make use of the exponential substitution ^ = exp(AT)\// which leads to an 
eigenvalue problem with transversality conditions 



M(Q)^Q = A^Q, Q.bQ = = QvQ; 
to make it short we can further omit the index Q 

(M-AU)^ = 0, Q-b = = Q-v. 



(77) 



(78) 



Should we substitute the incompressibiUty and transversaUty conditions 

avVy + QyVv , QA + Qyb, 

^z — :^ — t)z = — 1 (79) 

Qz Qz 

in the secular equation, we would end up with an overdetermined set of equations. To avoid it, we omit the equations 
which initially have Afe- and Xv^ terms, i.e., the 3-rd and the 6-th rows in the secular equation. In such a way we derive 
a secular equation for a reduced matrix 



(M-A]l)v/ = 0, = 



This secular equation 



by 



M : 





-Qa 







1 -v,;,e' 





-Qa 





\ 



Qa 






Qa 



2;vi.v(ffl+l)-v[e- 
2;7y;7y(a)+l)-(2a) + l) 



'2n,ny(0-v!^Q^ j 



(80) 



(81) 



A(A;Q,v^,Vk) = det(M-All) =0 

has 4 eigenvalues and via a calculation of the eigenvectors, we can derive h- and v- according to the transversality 
conditions given by Eq. (79). 

For an ideal fluid both v,J, and Vj^ are equal to zero. Omitting the viscosity terms we have a relatively simple form 
for the secular equation 



A(A;Q,v4 



0) 



2nyn^X^ - 



{[(4-8«>2+4_4„2] ^2^ [(2-8«;)«2-4«2 + 2] (a^+iaDX" 

- lQ\nyn,X + 22^ („2 + l ) (o^ + g4 ^ Q. (82) 

As we pointed out these eigenvalues give only a WKB approximation for the dynamics of MHD variables )/(t). For 
the special case of Qy = 0, which corresponds to an axial-symmetric motion, with a rotation along the z-axis, the 
secular equation gives directly the growth rates of the linearized MHD equations. 

F4(A ; e, = 0, v; = 0, = 0) = + 2 [e^ + ( 1 + 2a)J («? + l)co^]X^ + 2QI («? + 1 ) co, + = 0. 

The most restricted case is for the wave-vectors parallel to the rotation axis Q = Qe- when Qa — Qz cos 9 

= 0, Qy = 0, = 0, = 0) = A^ + 2 [e^ + (1 + 2a)Ja)J A^ + {qI + 2(o^)qI = 0. 

This is perhaps the most cited bi-quadratic equation in the whole history of science because it describes the magne- 
torotational instability (MRI) discovered by Velikhov [4] in 1959. In the astrophysics, this equation was recognized 
and overexposed by many astrophy sical grants 30 years later, see equation Eq. ( 111 ) of Ref . [2] and historical remarks 
therein. If we consider the special case of pure shear O)^ = with Qy = this dispersion equation gives the usual 



(83) 



(84) 



Alfven waves 



+ 2^)2^0, 



leal, 



(85) 



i.e., the rotation destabilizes the Alfven waves. The polarization of the magnetic field and the velocity of the Alfven 
waves are along the shear flow. 

For pure axial magnetic field B ^Be^, i.e., a = (0, 0, 1), and Qa ~ Qz - The matrix reduction is then given by simply 
erasing the z-components and taking into account only the x- and y-projections of the equations of motions 





/ 








-Qz 









1 








-Qz 


MRI — 




Qz 














v 





Qz 


-(20)^ + 1) 







(86) 



The secular equation is equation (84) for the MRI with Qa = Qz- 

The projection method can be generahzed in the general case if we introduce 2 unit vectors perpendicular to the 
wave- vector Cq — Q/Q 



- e2 



|l) = ei 



ICzXCqI 



e2 xeQ 
|e2 X cqI 




(87) 



(88) 



and also the corresponding bra-vectors 



Q.Qz,QyQz,~Ql-Ql) ^gg^ 



(- e,,av,o) 

Ql + Ql ' 



(90) 

For the degenerated case of Qx = = Qy we can regularize by choosing = i and Qy = 0. Then the limit i — )• gives 
the regularizations |1) = e.v and |2) = Cy. For all matrices Mq.^ where a, P =b, v we calculate the matrix elements in 
the two-dimensional space 

(Mo,.p).., = (j\Ma.p\}'): Where j,j' = 1,2. (91) 
In such a way we obtain a reduced 4x4 matrix 



(92) 



whose eigenvectors are automatically perpendicular to Q, simply because we have used the orthogonal to Q space. 
As a rule the linearized analysis is made in Lagrangian, moving, wave-vector space 

d,K(T)=Usheai(K(T)), (93) 

with a time-dependent wave-vector 

Kx = K^fi - -K^v(t - To), Ky = const, K- = const (94) 

for each MHD wave. 

In these coordinates for linearized waves the substantial time derivative D^'^''" = d^ is reduced to a usual time 
derivative and the separation of variables gives a set of ordinary independent equations for every MHD wave 

d,^K(T) = M(K(T))^K(T), K(t)-vk(t)=0, K(t) • = 0. (95) 

In this linearized case it is possible to exclude and V;. In such a way we arrive at a simple-for-programming set of 4 
equations and 2 zero divergence conditions 

d,v/K(T) = M(K(T))vAK(T), b, = -{K,{T)b, + Kyby)/K„ v, = -{K,{x)v, + KyVy)/K,. (96) 

For small Ky one can apply the WKB approximation supposing exponential time dependence of the MHD variables 
^(t) exp(AT) and the wave amplitudes. In the WKB approximation the energy amplification between T = —°° and 
T = +00 is given by the eigenvalue A with the maximal real part 

G « exp (2 dTReAn,ax(K(T)) j . (97) 

For the case of MRI with nonzero B- the amplification factors are so giant that the linear analysis makes no sense 
because the nonlinear terms become rather important and we have a nonlinear saturation of the MRI. This saturation 
simulates strong turbulence for small wave-vectors, but definitely for large wave-vectors l/Tyj ^ 1 at T ^ 00 we have a 
wave type turbulence with a given frequency. 

We have to mention that the linearized case of pure shear is exactly integrable in terms of the Heun functions [1,5]. 
Investigating numerically this case with 0)^, = and B- = in his Ph.D. thesis [6] T. Hristov discovered in 1990 the 
amplification of slow magnetosonic waves (SMWs) by shear flows. Applied to the physics of accretion disks this 
amplification works even for purely azimuthal magnetic fields and gives a scenario for weak magnetic turbulence 
related to amplification of SMWs. We had to wait 30 years of incubation period, cf. [7], for the SMWs amplification 
to be recognized as an important issue for that astrophysical phenomenon. In the next section we will consider how to 
proceed with the solution to the MHD equations. 



6. ENERGY DENSITY AND POWER DENSITY 

Our first step is to calculate the energy of plane MHD waves with time-dependent amplitudes. Using that 

J e^-^d^X = (2;r)35(Q) (98) 

for the energy we obtain 

\ J (pVLe + ^Bwave) = pV^^^, EE ^(v^^ + 5^^), (99) 

i.e., the energy density is 

Analogously, with the help of the viscous stress tensor d/^, we express the volume density of the wave heating 

<2kir = / CT45,yr™d3x cy;,{d,vr' + d^vrnd'x = f / idkVr' + d^Vm'-d'x (IOD 



^ i:av,e'Q-^+£e^v,e'Q'M d^x=pAylv(i:e^v^. 

\Q Q' / Q 



Similarly for the Ohmic part of the energy dissipation rate we have 

?2 



(102) 



gwave ^ j . E ^ ^ (rOtBwave)' = 4^ E V X bQe'^'^ = pAVlv^ ^ 

The dissipation rate of a laminar shear flow is given according to Newton's formula 

Q&r = \j (^*^*"" + diVf^'n^d^ ^ ^a2(5,,,5,,, + 5u5k,y? = riA\ (103) 
Now we can calculate the total energy dissipation gtot — G^n" + 2kin™ + Qci""^' ^^e viscosity and the effective viscosity 

J7eff 

e shear <^ 
~ ^2 ' I7eff — P Veff — . (i04) 

In this way we can express the effective kinematic viscosity by the dimensionless Fourier components of the velocity 
and the magnetic field 

Veff(T) = Vk + VkEe'v^(T) + V,nEG'b^(T)- (105) 

Q Q 

For example, if we have static probability distribution functions for the velocity and the magnetic field, the enhance- 
ment factor of the effective viscosity is given by the time-averaged squares of the Fourier components for T ^ 1 

'/ Q Q 

this important parameter determines the work of the accretion discs as a machine for making of compact astrophysical 
objects. The most simple scenario is to have the solution to the static equations for the i-th iteration of I'q and to 



calculate the next (/ + 1 )-th iteration 



('■+1) 



-Qy -.^ = -v-n 6,. + 2 — ^ — Q + 2(0 



(i+1) (i+i)x , , (i+i)„ ; 



(,■+1) 



-(a • Q)bi5+'' - vlG2vg+^> + n^Q^ [v« ® v« • Q' + (Q' x b«) x b«_Q, 

Q' 



dQ, 



(107) 



(108) 



5^ = 



-Qy^pr. T = Dt ^dr + df, for the independent variables (T,^^,,^;), = -QyT- (109) 

c'G.r Gy 



For cold protoplanetary disks the Ohmic resistivity of weakly ionized gas is very high and the effective viscosity is 
dominated in Eq. (106) by the Vm/Vk term, in other words, the viscosity of the protoplanetary disks is created by Ohmic 
dissipation. Completely opposite is the situation for the hot almost completely-ionized hydrogen plasma in quasars. 
The Ohmic resistivity is negligible and the effective viscosity is created by the Fourier components of the MHD waves 

^Vq^ . Only for small wave-vectors the MHD turbulence remains strong turbulence. At large wave-vectors we have 

weak wave turbulence with wave-vectors going to infinity. In the next section we will consider the stability conditions 
which have to be checked. 



7. STABILITY 

The linear Lyapunov analysis which we outlined in Sec. 5 gives the idea what we have to do when we obtain the 



static solution w'q = (bg'', Vq'). In order to investigate the stability of this static solution we have to consider a small 



time-dependent deviation from this solution 1'q''(t) ~ (bg' (t), Vq^ (t)). In this case, neglecting the quadratic terms 
with respect to 'Pq ■* , we find that the nonlinear terms in the MHD equations are linear integral operators in the Q-space 



,(1)/ 



V^j'iQ, + (Q' X bg,') X b^^ig, + (Q' X b^jV) X b^^ig, 



JO) 



-QxLQ,(v;j,'xb-Q, + v-xb-Q, 



(0) 



(110) 



We obtain new terms in the eigenvalue problem which finally is reduced to the problem of obtaining the maximal 
eigenvalue of an integral equation in which the coefficients are solutions to the static MHD equations. Now let us 
analyze the perspectives. 



8. PERSPECTIVE 

The open question of the missing viscosity in accretion disks is a longstanding problem in physics. In this work we 
have arrived at a complete set of numerically solvable ordinary differential equations which would help in finding a 
reasonable explanation of the aforementioned problem. 
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